Accretion disk in Cyg X-l in the soft state 
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Abstract. We analyze the simultaneous ASCA/ RXTE observation of 
Cyg X-l in the soft state on 1996 May 30. We apply a multi-color disk 
model with a torque-free boundary condition in the pseudo-Newtonian 
potential, and point out the main discrepancies between the commonly- 
used simplified 'diskbb' model and our one. We estimate the black hole 
mass and the accretion rate, and show that during this observation the 
disk probably did not reach the last stable orbit and the transition from 
the hard to the soft state was not yet completed. 



1. Introduction 

We study the cold disk spectrum of Cyg X-l in the soft state. The source was 
observed simultaneously by ASCA and RXTE on 1996 May 30. These data 
were previously analyzed by Dotani et al. (1997) and Cui et al. (1998). Dotani 
and his collaborators examined the ASCA data only. They used a multi-color 
disk blackbody model ('diskbb' in xspec, see e.g. Mitsuda et al. 1984) and 
found fcT max = 0.43 ± 0.01 keV and k dbh R in = 220+^km, where they give the 
correction factor for 'diskbb' model, fcdbb = 25/16 ~ 1.56. They also applied 
a general relativistic disk model ('grad' in XSPEC, Hanawa 1989), assuming the 
disk extends down to the last stable orbit in the Schwarzschild metric, and found 
a central object mass of M = \2^M.q. Cui et al. analysed simultaneous ASCA 
and RXTE observation. Their results concerning the disk parameters are similar 
to the above ones. 

2. The model 

The soft-state spectrum of Cyg X-l basically consists of two main components: 
the cold disk emission below ~3 keV and a power-law continuum extending 
into 7-rays (see Figure |]). We model the disk emission by the multi-color disk 
blackbody with the torque-free boundary condition in the pseudo-Newtonian 
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Figure 1. The simultaneous ASCA and RXTE observation of Cyg 
X-l in the soft state on 1996 May 30. The model consists of a non- 
thermal continuum, Compton reflection from the cold matter and a 
disk emission (see text for details of the disk model) . 



potential (see below). For the continuum we use a hybrid thermal/non-thermal 
model plus Compton reflection from the cold matter (Gierlihski et al. 1998). 
Our cold disk model is as follows. We consider a standard, flat, optically thick 
accretion disk around a non-rotating black hole of mass M. We use the following 
definitions: R g = GM/c 2 , R ms = 6GM/c 2 , r = R/R g , r = R/R ms and (3 = 
T co \/T c g. Subscript 'in' denotes the inner disk radius. We adopt the pseudo- 
Newtonian potential (Paczyhski & Wiita 1980), 
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which yields the color temperature distribution along the disk radius, 
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Figure 2. Temperature distribution along the disk radius, f = 
R/R ms , in a pseudo-Newtonian potential (lower curve) and in a sim- 
plified 'diskbb' model (upper curve). 

We create the model spectrum of the accretion disk integrating the observed 
flux, F v , numerically. Our model has three parameters: the inner disk radius, 
r; n = Rin/Rg = 6f; n , the maximum temperature, T max , and the normalization, 

from which the absolute value of the inner disk radius, i?i n , can be derived, 
assuming distance to the source, D, and the disk inclination angle, i. 

We point out that the above model differs from the commonly-used 'diskbb' 
model, in which the torque-free boundary condition is neglected and temperature 
distribution is simplified, 

T(f) = T f- 3 / 4 . (6) 

Discrepancy between the temperature distributions, as shown on Figure ^, is 
crucial when the disk extends down to the last stable orbit. Nevertheless, the 
integrated spectrum, F u , is similar, after some shift in parameters. We find 
that in the case of r ia = 6 our model spectrum for given R m and T max may 
be approximated by the 'diskbb' spectrum, for i?f n bb ~ 2.73i?; n and T^ b x w 
1.04iW 

The two models differ significantly in the accretion efficiency. The total 
disk luminosity in the pseudo-Newtonian disk (PN) model is 

L PN = ijtfc 2 , (7) 
lb 



giving the accretion efficiency, r^p^ = 1/16 = 0.0625, very close to the exact 
value in the Schwarzschild metric, rjs = 1 - (8/9) 1 / 2 ~ 0.0572 (e.g. Kato, Fukue 
& Mineshige 1998), whereas the 'diskbb' model yields 

L dhb = \mc 2 , (8) 

with very inaccurate efficiency, %iskbb = 1/4. Therefore, the accretion rate 
calculated from the 'diskbb' model is ~4 times less then the actual one. 

3. Disk extending down to the last stable orbit 

In all fits hereafter we assume D = 2.5 kpc, i = 45° and j3 = 1.7. 

First, we assume that the optically thick accretion disk extends all the way 
down to the last stable orbit in the Schwarzschild metric, r- m = 6. From the spec- 
tral fitting we obtain fcT max = 350 ± 10 eV and i?; n = 1481s 1 km, corresponding 

to a central object mass, M = 16.7^o;gM and an accretion rate, M ss 0.8 x 10 18 
g s _1 . 

The result obtained from the 'diskbb' model (Dotani et al. 1997), &dbb-Rin = 
220^ km, yields i?i n « 141 km, when fcdbb = 25/16 is applied. However, as we 
have shown in the previous section, fcdbb ~ 2.73, and this result should be rather 
i?i n 81 km, much less that we have found. The similar value can be computed 
for the paper by Cui et al. (1998). The discrepancy between R\ n derived from 
our model and these results comes from the different continuum model assumed. 
Dotani et al. (1997) used simply a broken power law, which fails to approximate 
properly the Comptonization spectrum at low energies. Cui et al. (1998) adopted 
a thermal Comptonization continuum, which at low energies should agree with 
our non-thermal model. However, the pure thermal model fails to reproduce 
the observed high-energy tail and cannot fit the data properly (Gierlihski et al. 
1998). Due to this limitation Cui and his collaborators obtained a very high 
temperature of seed photons for Comptonization, kT s ss 1 keV, which, in turn, 
affected the low-energy part of their model spectrum. 

4. Inner disk edge is farther? 

An interesting result can be obtained, when the relative inner disk radius, r- m = 
Rm/Rg-, is allowed to be a free fit parameter. The shape of the disk spectrum 
changes slightly with increasing rj n , and we find these changes large enough to 
constrain r- m by the data. A fit is improved (by A% 2 = 6.0/574) with the new 
best-fit value of r- m = 17tg 6 . The case of the disk extending down to the last 
stable orbit is outside the error bars here, suggesting that presumably Cyg X-l 
did not finish its transition to the soft state yet. This result is further supported 
by temporal properties of Cyg X-l. The evolution of power spectra, X-ray time 
lags and a coherence function (Cui et al. 1997) indicate that on 1996 May 30 the 
source was still on transition from the hard to the soft state, while in June the 
soft state was established. We note however, that unless all relativistic effects in 
the strong gravity are properly taken into account, the uncertainties of our disk 
emission and reflection models do not allow us to make any strong statements. 
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Figure 3. The solutions of the optically thick disk with a fraction of 
accretion power, /, dissipated in the corona, for m = 16.7 and a = 0.1, 
at a radius of the most significant instability R = 11.4R g (Svensson & 
Zdziarski 1995). Three curves from the left to the right correspond to 
/ = 0, 0.3 and 0.9. The horizontal line corresponds to the accretion 
rate of Cyg X-l in the soft state, rh « 0.5. 

The other disk parameters found in the rj n -free model are fcT max = 375j^eV 
and R in = 300l|fkm, yielding M = 12^M and M » 1.2 x 10 18 g s -1 . 

5. Stability of the disk 

To find out whether the cold disk on May 30 did reach the last stable orbit 
or not, we may try to investigate the stability of the disk. Let us consider a 
standard Shakura-Sunyaev disk with a fraction / of the gravitational energy 
dissipated in a corona. This model implies a critical accretion rate, m cr it, above 
which a radiation-pressure dominated region exists in the disk, which is thought 
to cause an instability. The critical rate is given by (Svensson & Zdziarski 1995) 

m crit = 0.37(am)- 1/8 (l-f)- 9/8 , (9) 

where m = M/Mq, a is the standard viscosity parameter and the relative 
accretion rate is defined as rh = Mc 2 /-LecM- 

Now we examine a case of the disk extending down to the last stable orbit, 
for which we found the mass, m = 16.7, and the accretion rate, rh = 0.5. From 
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the compactness ratio, lh/h = 0-36 (Gierlihski et al. 1998), we estimate / 0.3. 
For a = 0.1 the critical rate is m cr it ~ 0.52. The actual accretion rate locates the 
soft state of Cyg X-l at the end of a stable branch of the disk solution. Figure 
H shows the solutions of the optically thick accretion disk with the corona for / 
= 0, 0.3 and 0.9. 

On the other hand, if we drop the assumption of the disk extending down to 
the last stable orbit, we find m = 12 and m cr ;t ~ 0.54. The actual total accretion 
rate, m ~ 1, is well above the critical one. The inner part of the disk is thus 
unstable, which, in fact, is consistent with the truncation in the disk below r- m ~ 
17, we observe. Instead of a cold disk, an advective, optically thin, hot solution 
might exists in this region. The accretion below 17 R g dissipates about 0.3 of a 
total power released in the disk, but due to advection the luminosity of the inner 
flow might be considerably diminished. Therefore, its thermal emission would be 
difficult to disentangle from the observed thermal/non-thermal spectrum, and 
its existence hard to confirm or to rule-out. 
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